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ABSTRACT 



We discover that the very young very low-mass star ISO 143 (M5) is driving an outflow based on spectro-astrometry of forbidden 
[S II] emission lines at 6716 A and 6731 A observed in UVES/VLT spectra. This adds another object to the handful of brown dwarfs 
and very low-mass stars (M5-M8) for which an outflow has been confirmed and which show that the T Tauri phase continues at 
very low masses. We find the outflow of ISO 143 to be intrinsically asymmetric and the accretion disk to not obscure the outflow, as 
only the red outflow component is visible in the [S II] lines. ISO 143 is only the third T Tauri object showing a stronger red outflow 
component in spectro-astrometry, after RW Aur (G5) and ISO 217 (M6.25). We show here that, including ISO 143, two out of seven 
outflows confirmed in the very low-mass regime (M5-M8) are intrinsically asymmetric. We measure a spatial extension of the outflow 
in [S II] of up to 200-300 mas (about 30-50 AU) and velocities of up to 50-70 km s -1 . We furthermore detect line emission of ISO 143 
in Call (8498 A), OI (8446 A), Hel (7065 A), and weakly in [Fell] (7155 A). Based on a line profile analysis and decomposition 
we demonstrate that (i) the Ca II emission can be attributed to chromospheric activity, a variable wind, and the magnetospheric infall 
zone, (ii) the O I emission mainly to accretion-related processes but also a wind, and (iii) the He I emission to chromospheric or 
coronal activity. We estimate a mass outflow rate of ISO 143 of ~ 10~ 10 M Q yr~' and a mass accretion rate in the range of ~ 10" 8 to 
~ 10~ 9 M G yr~' . These values are consistent with those of other brown dwarfs and very low-mass stars. The derived M ml ,/M acc ratio 
of 1-20 % does not support previous findings of this number being very large (>40 %) for very low-mass objects. 

Key words, stars: low-mass - stars: pre-main sequence - circumstellar matter - stars: winds and outflows - stars: individual (ISO 217) 
techniques: high angular resolution 



1. Introduction 

While many details of the origin of brown dwarfs are still un- 
known, it was established in the last few years that during their 
early evolution brown dwarfs resemble higher mass T Tauri 
stars in main properties. Very young brown dwarfs (a few Myr) 
show surface activity, such as cool spots (e.g., Joergens et al. 
2003). There is evidence that brown dwarfs have disks from mid- 
infrared (mid-IR, e.g., Comeron et al. 2000; Jayawardhana et 
al. 2003; Luhman et al. 2008) and far-IR/submm excess emis- 
sion (Klein et al. 2003; Scholz et al. 2006; Harvey et al. 2012a, 
2012b). Many of these disks have been found to be actively ac- 
creting (e.g., Mohanty et al. 2005; Herczeg & Hillenbrand 2008; 
Bacciotti et al. 2011; Rigliaco et al. 2011), and several show 
signs of grain growth and crystallization (e.g., Apai et al. 2005; 
Pascucci et al. 2009). Furthermore, very young brown dwarfs 
rotate on average more slowly (e.g., Joergens & Guenther 2001; 
Joergens et al. 2003; Caballero et al. 2004) than their older coun- 
terparts (e.g., Bailer- Jones & Mundt 2001; Mohanty & Basri 
2003), which is indicative of a magnetic braking mechanism due 
to interaction with the disk. 

Jets and outflows are known to be a by-product of accretion 
in the star-formation process, and they have been observed for 
many classical T Tauri stars (CTTS, e.g., Ray et al. 2007 for a 
review). It has recently been shown that brown dwarfs and very 
low-mass stars (VLMS) are also able to drive T Tauri-like out- 



flows based on the detection of forbidden emission lines (FELs, 
e.g., Fernandez & Comeron 2001; Looper et al. 2010), spectro- 
astrometry of FELs (Whelan et al. 2005, 2007, 2009a, 2009b; 
Bacciotti et al. 2011; Joergens et al. 2012), and direct imaging 
in the CO J=2-l transition (Phan-Bao et al. 2008). To date, out- 
flows have been confirmed for six M5 to M8 type brown dwarfs 
and VLMS providing evidence that objects of a tenth of a so- 
lar mass to less than 30Mj up can launch powerful outflows: Par- 
Lup3-4(M5),pOph 102 (M5.50), ISO 217 (M6.25), LS-RCrA 1 
(M6.5), 2M1207 (M8), and ISO-Oph32 (M8). 

We present here the discovery of an outflow of the VLMS 
ISO 143 (M5) based on the spectro-astrometry of [S II] lines 
in high-resolution UVES/VLT spectra. This adds another ob- 
ject to the small sample of detected outflows for brown dwarfs 
and VLMS. We investigate the outflow properties of ISO 143. 
Furthermore, we show that ISO 143 exhibits several other 
activity-related emission lines and localize their formation sites. 

The paper is organized as follows: After a summary of the 
known properties of ISO 143 (Sect. [2]), the acquisition and anal- 
ysis of high-resolution UVES/VLT spectra is described (Sect.0. 
In Sect.|4l the results of the emission line study and the spectro- 
astrometric detection of the outflow is presented, followed by 
our conclusions on ISO 143 in Sect.|5] 



* Based on observations obtained at the Very Large Telescope of the 
European Southern Observatory at Paranal, Chile in program 080. C- 
0904(A) 082.C-0023(A+B) 



Spectral type pOph 102: K. Luhman, private communication. 
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2. The very low-mass star IS0 143 

ISO 1420 is a very low-mass star of spectral type M5 situ- 
ated in the Chal star-forming region (Luhman 2007). An esti- 
mate of its mass based on a comparison of effective tempera- 
ture and luminosity (7^=3 125 K, Lboi =0.088 L , Luhman 2007) 
with evolutionary models (Baraffe et al. 1998) yields a value of 
0.18 M . ISO 143 exhibits mid- and far-IR excess emission pro- 
viding evidence of a circumstellar disk. These detections were 
made with the InfraRed Array Camera (IRAC) on board the 
Spitzer satellite (3.6, 4.5, 5.8, 8.0 /im), the Multiband Imaging 
Photometer for Spitzer (MIPS, 24yum, Damjanov et al. 2007; 
Luhman et al. 2008), and the Photoconductor Array Camera 
and Spectrometer (PACS) of the Herschel mission (70, 160/itn, 
Harvey et al. 2012b). An IR spectrum of ISO 143 taken by the 
Sp/fzer/InfraRed Spectrograph (IRS, 5-36/mi) shows a weak 
10//m silicate feature (Manoj et al. 2011). ISO 143 was classi- 
fied as a class II source based on the spectral slope between 8 fim 
and 24 /an (Luhman et al. 2008; Manoj et al. 2011). The disk 
appears to be relatively flat and to have undergone substantial 
dust settling (Manoj et al. 201 1). Modeling of its spectral energy 
distribution (SED), including Herschel/PACS data at 70 //m and 
160 /mi (Harvey et al. 2012b; Y. Liu, pers. comm.), yields a total 
disk mass of 5xlO~ 6 M (assuming a standard gas-to-dust ratio 
of 100), a disk flaring index of less than 1.1, and an inclination 
less than 60°, with the most likely value being 15°-25°, i.e. close 
to face on, although inclination values 0°-50° are almost equally 
likely. ISO 143 shows Ha in emission with an equivalent width 
of 1 18 A (Luhman 2004), indicating that the disk is actively ac- 
creting. 



3. High-resolution spectra and spectro-astrometry 

We observed ISO 143 with the Ultraviolet and Visual Echelle 
Spectrograph (UVES, Dekker et al. 2000) attached to the VLT 
8.2 m KUEYEN telescope in three nights in 2008 and 2009 in 
the red optical wavelength regime at a spectral resolution A/AA 
of 40000. The spatial sampling was 0.182"/pixel. An observing 
log is given in TableQ] 

We used UVES spectra that were reduced, wavelength-, 
and flux-calibrated by means of the ESO UVES pipeline in 
order to study activity related emission lines (Ca II, O I, He I, 
[S II], [Fell]) regarding their line profile shapes (Sect. 14. It . 
Furthermore, a spectro-astrometric analysis of the detected for- 
bidden line emission in [S II] was performed based on a custom- 
made reduction and wavelength calibration procedure. After 
completing a standard CCD reduction of the raw data (bias and 
flatfield correction and cosmic-ray elimination), a row-by-row 
wavelength calibration of the two-dimensional spectra of indi- 
vidual echelle orders was done using the longslit package (fit- 
coords/transform) of IRAF. Finally, the sky was subtracted. We 
measured the spectro-astrometric signature in the resulting two- 
dimensional spectra by Gaussian fitting the spatial profile at each 
wavelength of both the FELs and the adjacent continuum, fol- 
lowing, e.g., Hirth et al. (1994a). The spatial offset in the FEL 
was then computed relative to the continuum. To increase the 
signal-to-noise ratio (S/N), which is the limiting factor in ap- 
plying spectro-astrometry to faint objects, the data were binned 
in the wavelength direction. To rule out spectro-astrometric arti- 
facts, which can be caused, for example, by an asymmetric PSF 
(Brannigan et al. 2006), we applied spectro-astrometry to a pho- 
tospheric absorption line (K I /17699) and demonstrated that the 

2 Simbad name: ISO-Chal 143 



spectro-astrometric offset in this line does not exceed 50 mas 
(apart from one outlier, see Fig.Q] top row). Further details on 
the spectro-astrometric analysis can be found in Joergens et al. 
(2012). 

Slit position angle. The UVES observations of ISO 143 were 
obtained within the framework of a high-resolution spectro- 
scopic study of young brown dwarfs and VLMS that is opti- 
mized for radial velocity work (e.g., Joergens 2008a). Therefore, 
the slit orientation was kept aligned with the direction of the at- 
mospheric dispersion, and, as a consequence, changed relative to 
the sky during the observations. Each 2x1700 sec exposure sam- 
ples a range of on-sky slit position angles (PAs) of about ±8°, as 
given in TableQ] Despite the slit PA varying during the exposure, 
the data nonetheless allowed the detection of outflows, as shown 
by our confirmation of the bipolar outflow of ISO 2 17 (Joergens 
et al. 2012) and by the spectro-astrometry of ISO 143 presented 
here. The spectra of ISO 143 were taken at mean slit PAs of 3°, 
10°, and 40°, which enabled us to investigate the observed out- 
flow extension as a function of the slit PA (see Sect. 14.21 ). 

Stellar rest velocity. A value of yo=16.9kms _1 was adopted 
for the stellar rest velocity of ISO 143, which is the mean 
value of the radial velocities measured for ISO 143 based on a 
cross correlation of many photospheric lines in the same spectra 
(Joergens 2012; cf. TableQ}. All velocities in the following are 
given relative to this value of Vo. 

4. Results 

4. 1 . Observed emission lines 

The observations of ISO 143 at high spectral resolution with 
UVES reveal line emission in Call, OI, He I, [SII], and [Fell], 
which are signatures of chromospheric activity, ongoing accre- 
tion, winds, and/or outflowing material. We studied their line 
profile shapes, measured their equivalent widths (EWs) and full- 
width-at-half-maximum values (FWHMs) in order to localize 
their formation sites in the system. For this purpose, we used the 
UVES spectra (Table[TJ after reduction, wavelength-, and flux- 
calibration by the ESO UVES pipeline. Furthermore, a polyno- 
mial fit to the continuum emission adjacent to each line was sub- 
tracted from the emission line regions. Several emission line pro- 
files of ISO 143 display complex shapes that obviously consist 
of more than one emission component. We decomposed these 
line profiles into its constituent components using the IRAF task 
ngaussfit. We fitted up to three Gaussian functions to the pro- 
files leaving the peak flux, velocity at the peak, and FWHM for 
each component as free parameters. We furthermore measured 
EWs for each line component. Table|2] lists the observed emis- 
sion lines, the determined peak radial velocities V, FWHMs, and 
EWs for each individual component. Details on the individual 
emission lines are presented in the following. 

Ca II IR emission. ISO 143 exhibits a broad, asymmetric, and 
variable Ca II IR emission line at 8498 A, which clearly consists 
of more than one component. Figure|2]displays this line at three 
different epochs in 2008 and 2009 in the top panels. We note 
that the two other lines of this IR triplet are not covered by the 
observations because of a gap in wavelength between the CCD 
chips of the two-armed spectrograph. A decomposition of the 
Call line profile into its constituent components demonstrates 
that it can be fitted by three Gaussian functions (Fig.|2]i: (i) a 
narrow emission component, which is centered close to zero ve- 
locity; (ii) a broader blue-shifted component; and (iii) a very 
broad red-shifted component. The narrow emission component 
(peak at velocities between -lkms -1 and 7kms~ 1 ; FWHM of 
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Fig. 1. Spectro-astrometric analysis of spectral lines of ISO 143 for three different observing times and slit PAs. Top row: photo- 
spheric line Kl A1699 as a test for artifacts, middle row: FEL of [S II] /16716, bottom row: FEL of [S II] /16731. For each case, we 
display the line profiles in the top panels, which are the average of 8 pixel rows in the spatial direction centered on the continuum 
after continuum subtraction, and the spectro-astrometric plots in the bottom panels, which show the spatial offset vs. radial velocity 
of the continuum (black asterisks) and of the continuum subtracted spectral line (red squares). Velocities are given relative to th§ 
stellar rest velocity V of ISO 143. 
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Table 1. Observing log of ISO 143. 



Date 


HJD 


Exptime 


Seeing 


Slit PA 


RV 






[sec] 


[arcsec] 


[dog] 


[kms" 1 ] 


2008 03 19 


2454544.69473 


2x1700 


0.67 


10.4 + 8.3 


16.315 ±0.082 


2009 01 15 


2454846.84705 


2x1700 


0.55 


3.1 ±8.2 


17.273 ±0.191 


2009 02 10 


2454872.87430 


2x1700 


0.62 


42.2 ±7.8 


16.986 ±0.398 



Notes. The listed entries are: observing date, heliocentric julian date 
at the middle of the exposure, exposure time, the averaged seeing cor- 
rected by airmass, slit position angle (PA), and radial velocity (RV, from 
Joergens, in prep.). The slit PA was not kept fixed during observations, 
and each 2x 1700 sec exposure samples a PA range of about ± 8 deg, as 
listed. 



9-21 kms -1 , EW of 0.1-0.7 A) is attributed to chromospheric 
activity giving its narrow profile centered close to the stars ve- 
locity. The variability of a few kms -1 seen in this line compo- 
nent (Table[2| but not in the velocity of the star (Table[TJ could be 
explained by rotational modulation caused by inhomogenously 
distributed activity regions in the chromosphere, such as plages 
(e.g., Huerta et al. 2008; Joergens 2008b; Prato et al. 2008 for 
T Tauri related activity causing velocity variability of a few 
km s -1 ; e.g., Berdyugina 2005 for Ca II emission produced in the 
chromosphere/ chromospheric plages). The broader blue-shifted 
component in Call (peak at velocities between -lkms -1 and 
-36kms -1 ; FWHM of 41-60 kms" 1 ; EW of 0.6-0.9 A) is sug- 
gested to be the signature of a wind that is expanding at a vari- 
able velocity. The very broad red-shifted emission component 
(peak at 26-63 kms" 1 ; FWHM of 94-165 kms" 1 ; EW of 0.6- 
4. 1 A), which is significantly variable in both the peak velocity 
(AV = 37 km s -1 ) and strength (A EW = 3.5 A), is likely formed 
in the magnetospheric infall zone of the system (cf. Muzerolle 
et al. 1998 for Call emission from infalling material in CTTS). 
The sum of the three individual Gaussians (Fig.|2l top panels, 
pink dashed line), is the overall fit to the Ca II profile. 

O I emission. ISO 143 has strong, broad, and variable emis- 
sion in the O I line at 8446 A (Fig. [2 bottom panels). This emis- 
sion is centered on red-shifted velocities in the two spectra 
from 2009 and on blue-shifted velocities in the spectrum from 
2008. We show that the line can be fitted by a single red-shifted 
Gaussian function for the two spectra in 2009 and that it can 
be decomposed into a blue- and a red-shifted component in the 
spectrum from 2008. The profile of the red component of this 
OI line has a very broad velocity distribution (FWHM of 106- 
121 kms -1 ) and a peak at velocities that are moderately red- 
shifted (9- 19 kms -1 ) but not as strong as seen in the red compo- 
nent of Ca II. The red component of O I is significantly variable 
in both the peak velocity (AV = 10 kms -1 ) and strength (A EW 
= 3.4 A). We attribute the emission in the red component mainly 
to disk accretion and, to a lesser extent, also to infalling material. 
The blue component that appears in this line in 2008 has approx- 
imately the same velocity (-42 kms -1 ) as the blue component in 
the Call line at the same epoch (-36 kms -1 ) suggesting a com- 
mon origin in a variable wind. 

He I emission. We detect a prominent He I emission line at 
7065 A (Fig.EBUl. The relatively narrow-band profile (FWHM 
of 26-35 km s _1 ) of this line being centered on V=0 and display- 
ing little variability in its strength (A EW < 0. 1 A), and the high 
temperature and density required to excite this transition indi- 
cates an origin in the chromosphere or the transition zone be- 
tween the chromosphere and the corona (cf. e.g. Sicilia-Aguilar 
et al. 2012). We note that a second weaker blended He I line at 



7065.7 A emerges in the red wing of this He I line at 7065.2 A 
(Fig.0. 

Forbidden [S II] emission. ISO 143 shows strong emission 
in forbidden lines of sulfur [S II] at 6716 A and 6731 A (Fig. Q} 
with a velocity at the line peak of 4 1 -50 km s - 1 , a FWHM of 27- 
33 km s -1 , and an EW of 0.2 A for [S II]^6716 and of 0.4-0.5 A 
for [S II]/i6731 (Table|2]i. The [S II] line profiles appear to vary 
little both in their strength (A EW < 0. 1 A) and peak velocity 
(AV <1 kms -1 , with the exception of the [SII]/i6716 profile on 
Feb. 2009; however, here the location of the peak was hampered 
by a low S/N). Spectro-astrometry of these lines (Sect. l4.2t re- 
veals that they are formed in an outflow at a distance from the 
central source of up to 30-50 AU. For both [S II] lines, only a 
red-shifted component is visible. This implies that the outflow 
of ISO 143 is intrinsically asymmetric and that there is either no 
or only little obscuration of the outflow seen in [S II] by the disk. 

Forbidden [Fe II] emission. We detect weak forbidden line 
emission of [Fell] at 7155 A and tentatively at 7172A (Fig.|4j). 
The peak of the [Fe II] Al 155 line seems to be around zero or 
low blue-shifted velocities; however, the poor S/N of this line 
hampers further quantitative analysis. If confirmed, the observa- 
tion of solely blue-shifted forbidden [Fe II] emission of a very 
low-mass outflow system that has a stronger red outflow compo- 
nent in [S II] would resemble the case of the [Fe II] emission of 
the brown dwarf IS0217 (Joergens et al. 2012). 




-200 -150 -100 -50 50 100 150 200 



V (km/s) 

Fig. 3. HeU7065 emission line in UVES spectra of ISO 143 
at different observing times in 2008 and 2009 after continuum 
subtraction. The flux is given in arbitrary units. This He I line at 
7065.2 A is blended with another, weaker He I line at 7065.7 A, 
which shows up in the red line wing at a velocity of about 
+30kms -1 . 
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Table 2. Observed emission lines of ISO 143. 



line 


slit PA 


date 


blue component 




green component 


red component 










V 


FWHM 




V 


FWHM 


EW 


V 


FWHM 


hw 




[deg] 




Tkm/sl 


Tkm/sl 

1 J\_lll/ j J 


[A] 


Tkm/sl 


Tkm/sl 

1 J\_lll/ j J 


rAi 
t^j 


Tkm/sl 


Tkm/sl 

1 J\_lll/ j J 


[A] 


rs TT1 




?nno ni is 

z.\j\jy ui lj 














40.8 


27.3 




(Kl 1 (S 44 A 1 


1 4 


?oo8 o^ i o 














41.9 


28.6 






42.2 


2009 02 10 














49.8 


30.8 


-0.2 


rs TTl 


3 1 


900Q 01 IS 














44.2 


31.7 


-0 4 


(fillet 8? A^ 


i n 4 


7008 0? 1 Q 

iwO UJ 17 














44.5 


26.7 


S 

U.J 




42.2 


2009 02 10 














44.6 


32.5 


-0.4 


Hp T 
LLC L 


1 1 
j . i 


9000 01 IS 








-0.5 


34.5 


-0.4 








f'706S 19 


10.4 


7008 03 1 9 








-2.3 


26.2 


-0.5 










42.2 


2009 02 10 








2.4 


28.1 


-0.5 








OI 


3.1 


2009 01 15 














8.8 


120.2 


-3.7 


(8446.36 A) 


10.4 


2008 03 19 


-42.1 


61.3 


-1.5 








13.5 


106.1 


-2.1 




42.2 


2009 02 10 














18.8 


121.2 


-5.5 


Can 


3.1 


2009 01 15 


-1.0 


60.0 


-0.9 


6.6 


8.5 


-0.1 


42.0 


164.8 


-1.7 


(8498.02 A) 


10.4 


2008 03 19 


-36.3 


51.1 


-0.6 


-1.3 


20.9 


-0.5 


62.8 


163.6 


-0.6 




42.2 


2009 02 10 


-18.4 


41.4 


-0.6 


5.9 


20.2 


-0.7 


25.7 


94.1 


-4.1 



Notes. The listed entries are: line with laboratory wavelength, slit PA, 
observing date, and for each individual component the velocity at the 
line peak, the FWHM, and EW. The laboratory wavelengths are taken 
from the NIST database (http://physics.nist.gov/asd3, Ralchenko et al. 
2011). The epochs are listed in the order of increasing slit PA, as in 
Table[3] for clarity. 





Call 8498, PA=3.1 deg, 2009 01 15 




Call 8498, PA-10.4 deg, 2008 03 19 




Call 8498, PA=42.2 deg, 2009 02 10 


1.4 

1.2 




1.2 




3.0 




1.0 




1.0 
0.8 




2.5 
2.0 




U_ 0.6 




0.6 
0.4 




1.5 




0.4 
0.2 




0.2 




1.0 
0.5 




0.0 




0.0 






YiJ^^^.J X L.. J^iSw^ 


-0.2 




-0.2 




0.0 





-200 -150 -100 -50 50 100 150 200 

V (km/5) 



OI 8446, PA^3.1 deg, 2009 01 15 



-200 -150 -100 -50 50 100 150 200 

V (km/s) 



OI 8446, PA^10.4 deg, 2008 03 19 






If \ HA 




// V" IV 









-200 -150 -100 -50 50 100 150 200 

V (km/s) 



OI 8446, PA^42.2 deg, 2009 02 10 




V (km/s) 



V (km/s) 



V (km/s) 



Fig. 2. Emission lines of Ca II /18498 (top row) and O I /18446 (bottom row) in UVES spectra of ISO 143 at different observing times 
in Mar. 2008, Jan. 2009, and Feb. 2009 in the order of increasing slit PA after continuum subtraction. The flux is given in arbitrary 
units. Gaussian fitting shows that the Call line can be decomposed into three components: (i) a narrow one centered close to zero 
velocity (probably chromospheric, green dotted line), (ii) a broader blue-shifted one (possibly produced in a wind, blue dotted line), 
and (iii) a very broad red-shifted one (possibly produced in magnetospheric infall, red dotted line). The pink dashed line is the sum 
of the Gaussian functions. The O I line can be fitted by a single red-shifted component in the two spectra from 2009 and by a blue- 
and a red-shifted component in the spectrum from 2008 (the pink dashed line denotes here the sum of the two Gaussian functions). 
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Table 3. Outflow parameters from forbidden [S II] lines and 
mass accretion rate from the OI 8446 line. 



Slit PA 


Date 


Line 


Offset 


v r 






Man, 


Lg446 


Mace 


[deg] 






[mas] 


[kms- 1 ] 


[cm" 3 ] 


Lb 


[M yr'] 


L Q 


[M yr'] 


3.1 


2009 01 15 


[SII] ,16716 
[SII] ,16731 


172 (±08) 
64 (±14) 


65 
54 


5600 


6.8xl0" 7 


3.1x10-'° 


3.6xl0~ 5 


1.8xl0" s 


10.4 


2008 03 19 


[SII] ,16716 
[SII] /16731 


138 (±17) 
95 (±16) 


48 

53 


> 2 x 10 4 


8.5xl0- 7 


8.3x10-" 


2.1xl0- 5 


8.0x10-" 


42.2 


2009 02 10 


[SII] ,167 16 
[SII] ,16731 


309 (±33) 
139 (±28) 


67 
52 


5600 


6.8xl0" 7 


3.0x10-'° 


5.4xl0- 5 


3.0X10 8 



Notes. The listed entries are: slit position angle (PA), observing date, 
detected FELs, spatial offsets of FELs in milli-arcsec measured by 
spectro-astrometry, corresponding radial velocities V r of FELs, electron 
density n e derived from the ratio [S II] ,1673 1/[S II] ,16716, line luminos- 
ity of [S II] ,16731, the mass outflow rate derived from [S II] ,16731 for 
an assumed outflow inclination angle of 45°, the line luminosity (of 
the red component) of 1,18446, and the mass accretion rate based on 

^8446- 



4.2. Results of spectro-astrometry 

We have clearly discovered the spectro-astrometric signature of 
an outflow in both [S II] lines at all three epochs in our UVES 
spectra of ISO 143. FigureQ] shows in the middle and bottom 
rows the line profiles of both the [S II] AA67 16,6731 lines and 
their measured spatial offsets as a function of the velocity. The 
spatial offsets are displayed for both the adjacent continuum and 
the continuum subtracted FEL. The plotted errors in the spectro- 
astrometric plots are based on lcr errors in the Gaussian fit pa- 
rameters. Table|3] lists the maximum spatial offsets and corre- 
sponding velocities of both [S II] lines in the order of increasing 
slit PA. Offset errors in Table[3] are based on the standard devi- 
ation in the continuum points. For an overview and to constrain 
the outflow PA (see below), we plot the maximum spatial offsets 
as a function of the slit PA in Fig. [5] 

Our spectro-astrometric analysis of the detected FELs of 
[S II] demonstrated that they originate from spatially offset po- 
sitions of the continuum source by up to 200-300 mas (about 



TiO Hel TiO TiO [Fell] [Fell] 




7060 7080 7100 7120 7140 7160 7180 
wavelength [A] 



Fig.4. [Fell] ^7155,7172 emission lines in a mean UVES 
spectrum of ISO 143. Also visible in this portion of the spectrum 
are TiO absorption bands with band-heads at 7053 A, 7087 A, 
and 7 124 A, and the He I XI '065 emission line. 



30-50 AU at the distance of Chal) at a velocity of up to 50- 
70kms- 1 . We found the [SII] ,16716 emission to be spatially 
more extended than the [SII] ,16731 emission (Table[3| Fig. [5), 
which is consistent with the [S II] /16716 line tracing lower den- 
sities than the [S II] /16731 line. 

The spectra were taken at mean slit PAs of 3° + 8°, 10° + 
8°, and 42° + 8° (cf. Sect.©. As illustrated in Fig.|5] we see 
a tendency for the estimated outflow extension to increase with 
increasing slit PA. The data thus hint that by varying the slit PA 
from 3° to 42°, one approaches the actual outflow PA. Follow- 
up observations, ideally spectra taken at two orthogonal slit PAs, 
are required to further constrain the outflow PA of ISO 143, in 
particular since the S/N is moderate for the spectra at 42°. 

4.3. Electron density and mass outflow rate 

The line intensity ratio [S II]^6731/[SII]/16716 can be used to 
derive the gas electron density n e in the line-emitting region be- 
cause these two lines have nearly the same excitation energy and, 
therefore, their relative excitation rates only depend on the ratio 
of collisional strengths, i.e. on the density (Osterbrock & Ferland 
2006, p. 121 ff; Bacciotti et al. 1995). We measured an [S II] line 
ratio for ISO 143 of 2.0 in spectra from Jan. 2009 (slit PA=3°) 
and Feb. 2009 (42°, cf. Tabled indicating an electron density 
« e of about 5600 cm 3 . For the epoch Mar. 2008 (10°), the esti- 
mated [S II] line ratio indicates an n e value above the line critical 
density n c , thus only a lower limit n e > « c ~ 2 x 10 4 cnr 3 can 
be derived in this case. 

We proceeded to estimate the mass loss rate of ISO 143 fol- 
lowing Hartigan et al. (1995; cf. also Comeron et al. 2003). For 
this purpose we derived the line luminosity of [S II]/16731 from 
its measured EW (Tabled and the dereddened Ro magnitude: 

L(L Q ) = 6.71 x 10~ 5 d 2 (pc) EW(A) 10" a4R ° . (1) 

We used a distance d of 162.5 pc (Luhman 2007), an R-band 
magnitude of 18.02 mag, which is obtained in the WFI R- 
band filter (Lopez-Marti et al. 2004), and an extinction of 
Ar(ivf/) =2.97 mag. The latter was converted from Ay=1.13 
(Luhman 2007) by interpolating the reddening law of Mathis 
(1990) for 7^=5.0 (Luhman 2004). The line luminosity we esti- 
mated for [S II]A6731 in ISO 143 is 7-9xl0 7 L (Tabled). 
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Fig. 5. Spatial offsets of FEL of ISO 143 as function of PA for 
[S II] /16716 (top panel) and [S II] A6731 (bottom panel). 



The mass loss rate was estimated from the line luminosity 
L6731 for n e < n c by the following equation (Comeron et al. 
2003): 



M ou t {Mq yr ) = 4.51 x 10 y 1 1 + ^ ) x 



L 613 i(L )V ± (kms ') 
l x (cm) 



(2) 



where n c is the critical density and V± and Z ± are the estimated 
tangential flow speed and scale length of the emission. In the 
case of the [SII]/16731 line that is within the high density limit 
(Mar. 2008, 10°), the mass loss rate does not depend upon the 
electron density (Hartigan et al. 1995) and Eq.|2]becomes 



M o „,(M yr- 1 ) = 4.51xlO 9 x 



Z ± (cm) 



(3) 



For the scale length of the emitting region l ± , we adopted the 
size of the aperture (cf. e.g. Hartigan et al. 1995), i.e. the UVES 
slit width of the observations of 1". For the outflow tangential 
velocity V±, we exploited the kinematic properties determined 
from spectro-astrometry of the [SII]/16731 line (Tabled, which 
indicate a radial velocity V r of the outflow of ISO 143 of around 
50kms _1 . Although the outflow inclination is unknown in prin- 
ciple, we can roughly estimate the tangential velocity V± from 
the observed V r by assuming that the outflow direction is perpen- 
dicular to the plane of the accretion disk and using the disk incli- 
nation value derived from the SED fitting (Harvey et al. 2012b; 
Y. Liu, pers. comm.). The disk inclination is, unfortunately, not 
very well constrained, because of the degeneracy of the SED 
model; the best-fit value is close to face-on (15°-25°), but disk 
inclination values 0°-50° are almost equally likely (cf. Sectf2]). 
The visibility of the red outflow component indicates an at least 
not completely face-on viewing angle of the disk, and we there- 
fore adopt, somewhat arbitrarily, an outflow inclination angle of 
45° for calculating the mass outflow rate, i.e. V r -V ± . 



The estimated outflow rate for ISO 143 (0.8- 
3xl0~ 10 M yr -1 , TableQ is consistent with previously de- 
termined mass loss rates of brown dwarfs and VLMS (Whelan 
et al. 2009a; Bacciotti et al. 2011). We note that assuming a 
larger outflow inclination angle, say of 65° (corresponding to 
a disk inclination of 25°), would result in outflow rates that 
are a factor of two lower than the given values. We find the 
same outflow rate for the observations at a slit PA of 3° and 
42° indicating no significant change as a function of slit PA 
from these data. While we derive a smaller outflow rate for 
observations at 10°, this deviation in the outflow rate by a factor 
of about four can be traced back to a small difference in the EW 
value of 0. 1 A and is, therefore, not taken as sound proof of a 
variable mass loss rate. Interestingly, a lower mass loss rate in 
this epoch would correlate with a decreased mass accretion rate 
in this epoch (see Sect. l4.4l i. 

4.4. Mass accretion rate 

The mass accretion rate of ISO 143 was derived from the accre- 
tion luminosity, 



R, 

Rin 



1.25 



~gm7 



(4) 



(Gullbring et al. 1998), by the use of the 01,18446 emission 
line as an indirect accretion tracer. We can assume that this 
line is predominantly produced by accretion, possibly except 
for the blue component we find in the profile of Mar. 2008 
(cf. Sect. 14. U . We determined the line luminosity of OI ,18446 
(for the spectrum of Mar. 2008 only the line luminosity of the 
red component) by applying Eq.Q] and using the dereddened i- 
band magnitude. We adopted an /-band magnitude of 15.51 mag 
(Luhman 2004), which was obtained in the DENIS /-band filter, 
and an extinction of A^denis) -231 mag, which is again con- 
verted from A j- 1 . 1 3 (Luhman 2007) by interpolating the extinc- 
tion law of Mathis (1990). A comparison of the determined Lg446 
(see TableO with an empirical correlation of line and accretion 
luminosity for OI ,18446 (Herczeg & Hillenbrand 2008) allowed 
us to roughly estimate the accretion luminosity L acc . This accre- 
tion luminosity is converted into a mass accretion rate by means 
of Eq.|4]by adopting a model-dependent mass of 0.18M o and 
a radius determined from the Stefan-Boltzmann law of 1.01 R e 
(cf. Sect.EJ. 

We estimate a mass accretion rate for ISO 143 between 
8xl0~ 9 (Mar. 2008) and 2-3xl0~ 8 M G yr" 1 (Jan. and Feb. 2009, 
Table[3]). These mass accretion rates are close to what was found 
for other VLMS, although at the higher edge of this distribution 
(e.g., Fang et al. 2009; Herczeg et al. 2009), but given the large 
scatter in the empirical Lu„ e -L acc correlation they are consistent. 
We can also convert the EW of Ha measured for ISO 143 (1 18 A, 
Luhman 2004) into a mass accretion rate by using the empirical 
relation of Fang et al. (2009) and assuming that the Ha emis- 
sion is entirely related to accretion. Following this path, we find 
a mass accretion rate of 1 x 10~ 9 M yr -1 , which is about an order 
of magnitude less than the value derived using the O I line. 

With the caveat in mind that accretion rate determinations are 
affected by large uncertainties, we find indications of a variable 
accretion rate of ISO 143 based on the O I line of 0.6 dex. This is 
consistent with other studies monitoring accretion rates of CTTS 
and that include VLMS, e.g. that of Costigan et al. (2012) who 
find variability of 0.4-0.9 dex for a sample of G2-M5.75 stars (cf. 
also Scholz & Jayawardhana 2006). The observation of accretion 
in ISO 143 derived from the OI line being on a lower level in 
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Mar. 2008 compared to Jan. and Feb. 2009 is also in line with 
the observations of the red component of the Ca 11/18498 line, 
which we attribute to accretion (cf. Sect. 14. H and Fig. [2]). 

Tentative hints were found of a relatively high ratio of mass- 
outflow to mass-accretion rate for brown dwarfs and VLMS 
(>40 %, Comeron et al. 2003; Whelan et al. 2009; Bacciotti et al. 
201 1) compared to that of CTTS (about 1-10 %, e.g. Ray et al. 
2007). The M 0Ut IM acc ratio we derive for ISO 143 appears to be 
closer to the CTTS value, as it is about 1 % using the O I-based 
accretion rate and about 20 % using the Ha-based accretion rate. 

5. Conclusions 

We have discovered that the young very low-mass star ISO 143 
(M5) is driving an outflow based on spectro-astrometry of for- 
bidden [S II] lines in UVES/VLT spectra. This adds another out- 
flow discovery to the very low-mass regime (M5-M8), where 
only about a handful of outflows have been confirmed so far. We 
have demonstrated that the forbidden [S II] emission in ISO 143 
is formed spatially offset from the central source by up to 200- 
300 mas (about 30-50 AU at the distance of Chal) at a veloc- 
ity of up to 50-70kms _1 . We found this outflow to be intrinsi- 
cally asymmetric as only the red outflow component is visible 
in the [S II] lines. This indicates an asymmetry in the launched 
outflow itself and/or in the distribution of gas in the immediate 
surrounding. Furthermore, in both cases there appears to be no 
or only little obscuration by the circumstellar accretion disk de- 
spite a low or intermediate disk inclination angle (< 60°) and 
a relatively flat disk geometry (Harvey et al. 2012b; Y. Liu, 
pers. comm.). To our knowledge ISO 143 is only the third de- 
tected T Tauri object showing a stronger red outflow compo- 
nent in spectro-astrometry, following RW Aur (G5, Hirth et al. 
1994b) and ISO 217 (M6.25, Whelan et al. 2009a; Joergens et 
al. 2012), although in general asymmetric jets are common for 
CTTS (-50%, Hirth et al. 1994b, 1997). We have shown here 
that, including ISO 143, two out of seven outflows confirmed in 
the very low-mass regime (M5-M8) are intrinsically asymmet- 
ric. We estimated the mass outflow rate of ISO 143 based on the 
[S II] lines to 0.8-3xlO~ 10 M G yr _1 , which is consistent with pre- 
viously determined mass-loss rates of brown dwarfs and VLMS 
(Whelan et al. 2009a; Bacciotti et al. 201 1). 

We furthermore detected line emission of ISO 143 from sev- 
eral activity -related emission lines of Ca II, O I, He I, and [Fe II] . 
Based on a line profile analysis we showed that likely forma- 
tion sites of these emission lines include the chromosphere, the 
accretion disk, the infall zone, and a wind. A decomposition 
of the Ca II /18498 IR line revealed three emission components, 
which can be attributed to narrow-band chromospheric emission, 
broader blue-shifted emission of a variable wind, and very broad 
red-shifted component probably produced in the magnetospheric 
infall zone of the system. The detected broad O I /18446 emission 
appears to be mainly formed in accretion-related processes in the 
disk. In addition, we observed a blue-shifted signature in the O I 
line at one epoch that is tracing the wind signature seen in the 
Ca II line suggesting a common formation site in a variable wind. 
The origin of the He 1/17065 emission can be attributed to the 
chromosphere or the transition zone between the chromosphere 
and the corona based on its narrow profile that is centered on 
the stars' velocity and on the high excitation temperature of this 
line transition. We estimated the mass accretion rate based on the 
OL18446 line to 0.8-3xl0~ 8 and based on the Ha EW from the 
literature (Luhman 2004) to 1 x 1 0~ 9 M yr~ 1 , which is consistent 
with that of other VLMS (e.g., Fang et al. 2009; Herczeg et al. 
2009). The derived ratio of the mass-outflow to mass-accretion 



rate M 0Ut IM acc of ISO 143 of 1-20% indicates a value close to 
that of other CTTS (about 1-10 %, e.g. Ray et al. 2007) and ap- 
pears to not support previous tentative findings of a very high 
M ou ,/M acc ratio for brown dwarfs and VLMS (>40 %, Comeron 
et al. 2003; Whelan et al. 2009; Bacciotti et al. 201 1). 

We conclude that the very low-mass star ISO 143 exhibits 
the well known activity features seen for higher mass T Tauri 
stars, such as chromospheric activity, accretion, magnetospheric 
infall, winds, and an outflow providing further confirmation that 
the T Tauri phase continues at very low masses. 
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